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Elastic Energy, Fluctuations and Temperature for Granular Materials 
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We probe, using a model system, elastic and kinetic energies for sheared granular materials. For 
large enough P/E y (pressure/ Young's modulus) and P/pv 2 (P/kinetic energy density) elastic dom- 
inates kinetic energy, and energy fluctuations become primarily elastic in nature. This regime has 
likely been reached in recent experiments. We consider a generalization of the granular temperature, 
T g , with both kinetic and elastic terms and that changes smoothly from one regime to the other. 
This T g is roughly consistent with a temperature adapted from equilibrium statistical mechanics. 

PACS numbers: 45.70.-n,05.40.-a,64.70.pf 



We explore the role of elasticity in the energy and en- 
ergy fluctuations of sheared dense granular systems. For 
dilated gas-like granular states, energy fluctuations are 
frequently described in terms of a temperature, defined 
as the fluctuating part of the kinetic energy, Tfc = m < 
v 2 > /2. Here, v is the local random component of 
the velocity. This definition is predicated on assump- 
tions such as molecular chaos, absence of correlations, 
and short-lived collisions, that do not always apply. For 
dense systems, a very different concept, Edwards entropy, 
has been proposed jlj . This quantity is the logarithm of 
the number of jammed configurations consistent with all 
constraints on the system, and the Edwards temperature 
is T^ 1 — OSe/QV, where V is the system volume. 

Both of these pictures assume that minimal energy is 
stored in compressional modes of the particles. This as- 
sumption is valid when the pressure is small compared 
to the Young's modulus, E y . However, there are situa- 
tions when this need not be the case. The main goal of 
this letter is to analyze via discrete element simulations 
(DES) the storage of energy and energy fluctuations for 
dense granular material subject to shearing. 

In order to better establish a context, we estimate the 
relative importance of kinetic and elastic energy for a 
simplified system. We imagine that spherical particles 
with a typical velocity, v, are subject to an applied force, 
L, at each contact. For simplicity, assume at first just 
a pair of opposing contacts, so that there is an effective 
pressure, P = L/A, where A = ttR 2 , and R is the radius 
of a sphere. The elastic energy per contact (assuming 
a Hertz-Mindlin contact law) is e oc L 5 / 3 . The ratio of 
elastic to kinetic energy for a particle is Re — e/K — 
C(P/E y ) 2 / 3 (P/pv 2 ), where C is an 0(1) constant that 
depends on the Poisson ratio, E y is the Young's modulus, 
and p is the bulk density. By multiplying by the number 
of contacts per particle, it is possible to generalize this 
to a more realistic situation. When Re is small, elastic 
energy is irrelevant, and vice versa when it is large, it 
should be included in a description of the system. 

We now consider Re for two representative cases. For 



the experiments of 0, E y = 5MPa, typical P ~ 160Pa, 
and p = 1.2 g/cm 3 . Typical speeds ranged over 6 x 
10- 4 m/s <v< 2xl0" 2 m/s, so that 0.3 < R E < 4xl0 2 . 
For glass spheres @, E y is larger by roughly a factor of 
5000. Assuming densities p ~ 2 g/cm 3 , and pressures 
corresponding to the base of a column 10 cm high in a 
gravitational field, Re — 1 for velocities of a few mm/s. 

Hence for slowly sheared dense granular systems, there 
exist velocity regimes for which elastic energy is the dom- 
inant mode of energy storage. In such a setting, neither 
Tfc nor Te is likely to provide a good measure of the ran- 
dom nature of the system. 

In this context, we propose an extension of granular 
"temperature" that contains information on fluctuations 
of the elastic energy, and then compare this extension 
to a relation for temperature drawn from statistical me- 
chanics. (A somewhat similar approach of applying equi- 
librium statistical theory has been taken in recent works 
on foams and granular systems p|.) To carry out this 
exploration, we use DES of 2D particles that are subject 
to plane shear and (possibly) compression. 

The generalization of "temperature" that we consider 
is based on the classical idea that for a lattice of elas- 
tic particles, the average fluctuating energy/particle is 
3fcsT. Using this as a heuristic guide, we define a gen- 
eralized temperature that is roughly T g = m < v 2 > 
/2 + k < x 2 > /2, where v corresponds to the fluctuating 
part of the velocity, and x to the fluctuating part of the 
compression of a particle. Note that this definition pro- 
vides a simple bridge between the extremes of a gas-like 
state and a highly compressed slowly evolving state. 

In the simulation, particles are confined between two 
inpenetrable straight parallel boundaries, as sketched in 
Fig. The top boundary, which is 50 mean parti- 

cle diameters (d m ) long, moves at a steady speed, and 
induces shearing in the system. The boundary condi- 
tions in the shearing direction are periodic. This system 
avoids the nonuniformity that characterizes physical ex- 
periment, which typically exhibit shear bands Q. Al- 
though we could have chosen to use the even simpler 
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Lees-Edwards conditions, we have concentrated on the 
present model so as to explore the influence of bound- 
aries, typically present in physical experiments. 

These simulations closely follow the soft-disk/sphere 
model to describe various granular systems (see Q 
and references therein). Here, we concentrate on two- 
dimensional polydisperse disks in a zero-g environment. 
The walls are made of identical particles that are rigidly 
attached. Forces between the particles have a normal 
component given by F^v = [kfX — 7jvW2.(v£j • n)] where 
kf is a force constant, r£j = \ri t j\, r£j = i"; — Tj, 



dij are the diameters of 



n = Tij/r^j, d = (di + dj)/2 
the particles i and j, x = d — r^j is the compression 
v ij" — — Vj , fh is the reduced mass, and 7™ is the damp- 
ing constant related to the coefficient of restitution, e n . 
The parameters represent photoelastic disks 0; in par- 
ticular, 7„ corresponds to e n = 0.5. The tangential force 
is given by F s = sign{~vl el )min{^ s fh\v t rel \,v k \F c N \)s 
where i>* ei is the relative velocity in the tangential di- 
rection s, 7 S = 7n/2 and /Life is the coefficient of friction 
between the particles. The equations of motion are then 
integrated using a 4th order predictor-corrector method. 
Additional simulations (to be presented elsewhere) show 
that variation of the parameters or the force model mod- 
ify only details of the results. 

The simulations are performed with approximately 
2000 polydisperse particles, with a radius variability of 
10%. While polydispersity is important to avoid crys- 
tallization, the details of the size distribution are not: 
the results are very similar for different ranges of particle 
sizes, or a bidisperse distribution. Particles are initially 
placed on a lattice, given random velocities, and the sys- 
tem is then very slowly compressed to a desired volume 
fraction, v. The results that follow use t, the time it takes 
the shearing wall to travel once across the domain, as a 
time scale, and l/t (I = 50<i m ) as a velocity scale. 

The quantities below are calculated using space-time 
averaging. Thus, the system is divided into cells, and 
averaged quantities are calculated for each cell. In par- 
ticular, the kinetic temperature is defined by 



n = l 



(m(..') 2 ) + (m(»') J ) + j{m{d r Jf) 



where u, v are the components of particle velocity, d p 
is the diameter of a particle, m is its mass (m ~ dp), 
(3 = 1/4 for disks, and u> is the angular velocity. The 
primed averages are defined with zero mean, e.g, ivfv') = 
(uv) — (u)(v). 

The elastic energy is obtained by averaging per colli- 
sion, not per particle. The difference between the two 
is significant for dense granular systems considered here, 
since particles typically experience multiple collisions. If 
the compression of particle j due to the collision 
c, then our definition of the elastic energy in cell / is 
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FIG. 1: (a) The geometry of the simulations (only centers 
of the particles are shown). (b,c,) Elastic energy and kinetic 
temperature, scaled by v 2 , as functions of the distance from 
the shearing wall, y. The shearing wall is to the right, and 
the shearing is in the out-of-plane direction (v = 0.1). 



where n; is the number of particles in cell I at a given 
time, and n; is the average number of particles during the 
period of Nt 3> 1 time steps (in practice, the averaging 
time scale is sufficiently short so that to a high degree 
of accuracy ni — hi). Definition Q ignores the energy 
used to overcome friction, and its form clearly depends 
on the form of the normal force, e.g., for a 3D systems, 
for which a Hertzian interaction law (Fjv ~ x 3 ^ 2 ) is more 
appropriate, E ei i is of different form as well. 

Figure n shows the elastic energy and kinetic temper- 
ature (scaled by the square of the shearing velocity v, 
and by the average mass of a particle) vs. distance, y, 
from the shearing wall. For these simulations, v is con- 
tinuously increased by (slow) compression, v increases 
from 65% (bottom) to 90% (top) (note that in 2D, ran- 
dom close packing and cubic close packing correspond to 
about 85% and 2n/y/3 ~ 90%, respectively). Clearly, 
as v is increased, there is a transition region (about 
v c = 80%) where the energy stored in the internal degrees 
of freedom (elastic energy) becomes more relevant than 
the kinetic energy (see also Fig. |2J . The y-dependence of 
the results is rather weak and becomes even weaker for 
higher i/s. Hence, hereafter, we ignore the y-dependence 
and use system averages of locally computed quantities. 

As v increases, the energy is mainly elastic, and Tj. 
loses its relevance. In order to have a quantity that might 
play the same role as in a dense granular system we 
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propose a generalized granular 'temperature' by 

Tg = Xfc + T e , 



(2) 



as a sum of Tfc, and the 'elastic' part T e . T e is defined as 
the mean fluctuation of elastic energy, in a manner simi- 
lar to Tk , which is the mean fluctuation of kinetic energy. 
This definition follows the classical statistical mechanics 
result where the mean fluctuations in the combined elas- 
tic and kinetic energy of an oscillator are proportional 
to the temperature. However, unlike this classical case, 
there is no reason to expect equipartition between elastic 
and kinetic modes. Rather, the ratio T e /(T e + varies 
from in the dilute limit, to f in the dense limit. 

The definition of T e requires some care due to multi- 
ple collisions. One simple and natural definition that is 
consistent with the classical statistical definition of tem- 
perature for an oscillator is as follows. We first define the 
average elastic energy per particle in cell I as 



(E e ,i) = ^n c (xi) 2 = ^-n c 



n in 



N t mn. 



fe=l j=l c=l 



(3) 

where (xi) is the average compression per collision, and 
n c is the average number of collisions per particle. Then, 



T e .l = 



(Sx 2 ) = -^-n c ({x j:C - ( Xl )) 2 ) = E Etl - (E e ,i) 



(4) 

where the last equality easily follows using and ©. 

Figure shows T e , and (E e ) vs. time, for four 
different v's. Unlike the results of Fig. ^ these results 
are obtained after shearing for long times at fixed u's. 
For higher v's, clearly T e ^> Tf., Interestingly, there are 
rather large fluctuations of the results with v just above 
v c rj 80%, corresponding to the regime where elastic en- 
ergy becomes predominant. Also, experiments have indi- 
cated a phase transition for comparable densities which 
may be related (see also 0). Generally, one might 
expect both glassy and/or jamming phenomena to dom- 
inate this regime, a point that we will explore elsewhere. 

We further interpret T g by comparing it to an alterna- 
tive definition from statistical mechanics. We consider, 
among various possibilities, the standard relation Q 



dU 
dT 



6U 2 

J>2 



(5) 



where U is the total energy in the (usually conservative) 
system, T is the temperature, and SU 2 = (U 2 ) — (U) 2 . 
We use U = E-^ in + (E e ), and ask whether it makes 
sense to define T = T g = T k + T e . To check this 
idea, we now think of (J5J as a defining equation for 
T m = VslP/{dU/dT g ) ( here T m stands for 'model tem- 
perature'). The agreement between T m and T g will pro- 
vide some indication of the utility of the definition for 
T g . Note that we should not expect perfect or possibly 
even any agreement, since we are considering a strongly 




FIG. 2: Elastic energy, kinetic and elastic temperatures 
various u's, scaled by v 2 . 
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dissipative system far from equilibrium. We emphasize 
that there are no fitting parameters and that the results 
that follow are obtained directly from simulations. 
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FIG. 3: Generalized temperature versus average energy for 
two different */s (increased U corresponds to increased shear- 
ing velocity v). The lines are the least square fits to the data. 
In the inset of b), U (solid line, filled circles) and T g (broken 
line, squares) are plotted versus ln(v) for slow shearing. 
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FIG. 4: Generalized, kinetic, and 'model' temperature (see 
the text and Eq. JSJ) for two volume fractions. Error bars 
(resulting from statistical uncertainty of the results, as well 
as from the uncertainty introduced by the least square fits 
from Fig. [3J are approximately represented by the size of the 
symbols. The insets show kinetic energy (filled circles) and 
average elastic energy (squares). 

Figure [3] shows typical data used to determine c v — 
dU /dT g needed for T m . The range of C/'s shown in Fig. El 
corresponds to 0.0001 < v < 40. For v = 80%, U ~ T g 
over almost 10 decades. For v — 86%, U ~ T g , where 
a ~ 2.0. However, in this case, there is a deviation from 
the power law fit for slow shearing, and we are limited to 
a smaller range of U's, since most of the energy is stored 
in the system as elastic energy. In order to obtain better 
data for U(T g ) for slow shearing at higher densities, we 
determine U and T g as functions of v, see Fig. |3Jd. We 
then find to a good approximation that both U and T g , 
are proportional to ln(v). This logarithmic dependence 
for slow shearing is, to the best of our knowledge, the 
first computational confirmation of recent experimental 
results |9j, and will be presented in more detail elsewhere. 
For our purposes here, it is sufficient to extract U(T g ). 

The dependence of U (T g ) is interesting. We note that 
U ~ T g for v = 80% is similar to a recent simulation 
of sheared foams 0. However, the result U ~ T g for 
v = 86% and for not too slow shearing (U > 1, see 
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Fig. |3Jd) is striking and deserves some comment. Al- 
though the increased role of elastic energy is important 
here, this does not appear to be the whole story, as 
seen in the context of Fig. Here, we show that even 
for this high v, for fast enough shearing the kinetic en- 
ergy is still dominant. Our preliminary interpretation is 
that increased volume fraction contributes significantly 
to decreased mobility (jamming) of the granular system, 
therefore reducing the increase of T g with U. 

Figure ^ contains a summary of the various types of 
temperatures considered in this study for two z/s. As 
already noted, for v — 80%, most of the energy is still 
kinetic, while for v — 86%, elastic energy is essential. 
For v — 80%, Tk is dominant, and also it satisfies the 
model J5J| since Tk T m . However, for v = 86% and for 
slow shearing, Tk is smaller by 4 orders of magnitude than 
T g and T m . Thus, T/. cannot be used to even approxi- 
mately describe a dense slowly sheared granular system. 
This difference decreases for higher shearing rates, but 
only at very high shearing is there good agreement. 

The agreement between T g and T m is not perfect, al- 
though prefect agreement is not to be expected. These 
studies demonstrate the clear need to incorporate elas- 
tic energy and elastic fluctuations, and that T g has util- 
ity as a generalized granular temperature. Nevertheless, 
there remain many open questions regarding the extent 
to which the various temperatures serve similar functions 
to their molecular counterpart. We will present detailed 
results addressing this type of question elsewhere. Here 
we note that at least qualitatively one can show that 
these temperatures can be used in the context of ther- 
mal conduction, i.e. the there is a flow of heat from hot to 
cold regions. We also note that the distributions of both 
kinetic and elastic energy are strongly non-Gaussian for 
a significant range of the parameters analyzed here. 

There are other possible tests of the proposed gener- 
alization, as recently discussed on the context of sheared 
foams IH- The ultimate test will be to perform physical 
experiments where the validity of the proposed concepts 
can be verified. In the theoretical direction, it will be 
of interest to relate the generalized granular temperature 
proposed here to the one resulting from Edwards-entropy 
based approach. 
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